Lande g-tensor in semiconductor nanostructures 
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Understanding the electronic structure of semiconductor nanostructures is not complete without 
a detailed description of their corresponding spin-related properties. Here we explore the response 
of the shell structure of InAs self- assembled quantum dots to magnetic fields oriented in several 
directions, allowing the mapping of the g-tensor modulus for the s and p shells. We found that the 
g-tensors for the s and p shells show a very different behavior. The s-state in being more localized 
allows the probing of the confining potential details by sweeping the magnetic field orientation from 
the growth direction towards the in-plane direction. As for the p-state, we found that the g-tensor 
modulus is closer to that of the surrounding GaAs, consistent with a larger delocalization. These 
results reveal further details of the confining potentials of self-assembled quantum dots that have not 
yet been probed, in addition to the assessment of the g-tensor, which is of fundamental importance 
for the implementation of spin related applications. 



Electronic magnetism in semiconductor nanostruc- 
tures is one of the important properties to be harnessed 
in spintronic devices [l| as well as in prototypical sys- 
tems for quantum information processing '2|. In order to 
understand and separate the effects of quantum confine- 
ment and band structure, including spin-orbit coupling, 
strain and non-parabolic effects, the response of the elec- 
tronic spin on an applied static magnetic field can pro- 
vide an improved picture of the overall quantum system. 
The electronic g tensor, which describes the symmetries 
and magnetic response of the unpaired electron system, 
is thus a very important tool to assess and investigate 
these fundamental aspects of spin electronics in nanos- 
tructures. 

For conduction electrons in bulk semiconductor crys- 
tals, the g-factor can be determined accurately by second- 
order k • p theory using Roth's equation[3, Q, and con- 
firmed by experiment [5]. For unpaired electrons bound 
to donors, g-tensor differences from the free atom value 
of 2 for the ground state will reveal the dependence on 
the crystal field and spin-orbit couphng[6|. The sym- 
metries of defects and chemical environment can be also 
revealed by mapping the g-tensor jTj. In addition to that, 
the anisotropic part of g influences spin-lattice relaxation 
and is important for spin- related applications [^, |9| . For 
the case of quantum wells jlfl] and wires [llj the g-tensor 
will be affected by quantum confinement, strain and com- 
position fluctuations. 

Experimental investigations of g factors and g te nsors 
have been reported for metallic nanoparticles 1211 and 
lithographically defined quantum dots (QDs) [12 



For metallic nanoparticles, the difficulties arise in find- 
ing the symmetry axis, which can be determined from 
the g-tensor mapping. In addition, because the electron 
mean free path is smaller than the particle size, angu- 
lar momentum may not be a good quantum number. 
For lithographically defined quantum dots, the Zeeman 
splitting and the orbital splitting have comparable energy 



scales, thus preventing the evaluation of the out-of-plane 
g-tensor component. For self-assembled quantum dots, a 
number of experiments have demonstrated striking sim- 
ilarities with the atomic behavior, such as Hund's rules 
and the Aufb au p rinciple in determining the shell filling 
for electrons [l5|, [la]- The charging, Zeeman splitting 
and single-particle energies are all different for this case, 
which allows them to distinguished for excitons and elec- 
trons 171, lla llSl, l20| . More recently, calculations were 
carried out displaying the relationship between the g ten- 
sor and the electronic structure for quantum dots 21 1 . 

In this Letter we explore the shell structure dependent 
spin properties of electrons trapped in InAs quantum 
dots (QDs). By evaluating the electron addition ener- 
gies inferred from magneto-capacitance data, we present 
an experimental account on the g-tensor modulus for the 
s and p states which were mapped out according to the 
crystallographic directions of highest symmetry. 

InAs QDs were grown by molecular beam epitaxy and 
capped with thin InGaAs strain reducing layers, as de- 
scribed elsewhere [l9|, |20|. These structures were em- 
bedded in capacitance structures that were subsequently 
defined by conventional photolithography. The area of 
the devices was 4 x 10~^mm^, hence encompassing an 
ensemble of about 10® QDs per device. Magnetocapaci- 
tance experiments were carried out at 2.7K for magnetic 
field intensities ranging from to 15T. Field sweeps were 
performed at 15° intervals covering at least 180° by tilt- 
ing the sample with a goniometer. 

Figure 1 (a) shows the second derivative of magneto- 
capacitance spectra taken for field sweeps along the [001] 
and [110] directions (polar scan). The energy scale de- 
rived from the applied bias and voltage-dependent lever 
arm is translated into the chemical potential within the 
QDs referenced to the GaAs conduction-band edge. The 
lever arm was calculated taking into account depletion 
effects in the back contact [20J, thus allowing the deter- 
mination of a precise energy scale. The gray scale map 



is a representation of the density of states (DOS), where 
the OD (up to - -80meV) and 2D (above - -80meV) 
levels associated with the QDs and wetting layer can be 
easily identified. For in-plane field all orbital effects are 
minimized given the pancake geometry of the quantum 
dots; also negligible are the effects of the magnetic field 
on the wetting layer, which exhibits Landau-level fillings 
for magnetic fields perpendicular to the sample surface. 




15 

Energy(meV) Magnetic Field (T) 

-222-217-212-207-202-197 „ 

> 

E 




(c) 


I^ 


^ 


"Tj ^ I OS-shell 
np-she!l 



-1.32 -1.27 -1.22 -1.17 

Voltage(V) 



3.0 6.0 9.0 12.0 15,C 
IVlagnetic Fieid (T) 



FIG. 1: (a) Measured data (derivative) for polar scan, i.e., 
magnetic field sweeping from parallel to the [001] direction 
towards the [110] direction. The grayscale is keyed to the 
second derivative of the capacitance spectra, with light colors 
indicating a higher density of states. Each of the 13 frames 
corresponds to B sweeps from to 15T, taken at 15° intervals 
with 9 the angle between B and [001]; The s, p and d states, 
as well as the wetting layer 2D levels are indicated as well; 
(b) Different capacitance spectra showing the filling of the 
s shell with 1 and 2 electrons, dots correspond to B = OT 
and triangle to |B| = 15T, and 9 — ■k/12. The solid lines 
correspond to two gaussian fits, (c) Fits to equation [l] for s 
and p shell electron addition energies for 9 = 7r/2. 

From the electron addition energies [16|, \l^ |20|, or- 
bital, electrostatic and Zeeman contributions can be sep- 
arated. Figure 1 (b) shows the shift of the capacitance 
peaks for the sequential charging of the s-level on the 
applied magnetic field for = 7r/12. The effects of 
Coulomb charging, diamagnetic shifts and Zeeman split- 
ting can be seen in these data. One can separate the 
orbital effect contribution from the others by subtract- 
ing the peak positions, as the diamagnetic contribution 
is the same for both. For the s shell, using a Fock- 
Darwin (FD) formalism 15|, |20| one has for the loading 
of the first and second electron for T=0 and B || [001]: 



where E^ is the confinement along the growth direction, 
51 = ^/Lu^+lJ[/4:, with tOc as the cyclotron frequency, 
E^ the Coulomb charging energy at zero magnetic field 
for the s shell electrons, and jg^^j/SB the Zeeman split- 
ting. For the s shell, one finds hojQ = 37.8±0.2meV and 
Ez ^ —280meV . Detailed modeling can provide an even 
better description [Ifi], however it is not necessary to 
capture the essential features. Electrostatic effects can 
be calculated given the single particle energies within 
the FD framework |15|]. For the s and p shells we find 
i?ss=17.2meV and £'pp=13.4meV. Figure 1 (c) shows the 
dependence of the addition energies for the s-s and p-p 
configurations on the applied magnetic field for 6 = 7r/2. 
We find that an agreement within 5%-15% can be found 
from the calculated and measured addition spectra. 

Two effects must be considered when analyzing these 
data - i) wavefunction compression and its effect on the 
charging energies and ii) temperature. As the mag- 
netic field is raised, the wavefunction is compressed 
which increases the Coulomb charging energies. Un- 
der the FD formalism, this effect can be calculated as 
El,{B) = ^^(0)(l-^a;2/4^2^i/4^ ^^^^^ ^i^^^^ j^ ^^^ 

Coulomb charging energy at zero magnetic field for the 
j-shell. This effect takes place for all directions of the 
applied magnetic field, but to a lesser extent for the in- 
planc configuration due to a stronger confinement along 
the growth direction. As a zeroth-order approximation, 
we assume this effect to be the same for all configurations. 
The consequence of this assumption is to underestimate 
the Zeeman contribution for in-plane magnetic fields. 

As stated above, the experiments were carried out with 
QDs ensembles. Thus, an accurate description of this sys- 
tem at finite temperatures requires usage of a magnetiza- 
tion model for a system of n non-interacting spins, which 
can be carried out by calculating the partition function 
for the system. For the current analysis, we take into 
account the temperature-dependent spin contribution on 
the addition energy spectra. This yields a direct relation- 
ship between the addition energy A/i, the B-dependent 
Coulomb charging, the Zeeman splitting and the temper- 
ature: 
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Figure 1 (c) shows A/i as a function of the magnetic 
field and the corresponding fit to equation [T] for the s 
shell dX 6 — 7r/2 [i.e., in-plane magnetic field). For the 
p shell, the same description applies, and a similar re- 
lation can be derived. However, the observed addition 
energies and the dependence on B are smaller (fig. Ic), 
and it becomes difficult to implement a reliable fit be- 
cause of two additional factors: i) a smaller Coulomb 
charging energy, and ii) broadening of the capacitance 
peaks due to the Coulomb disorder. Coulomb charging 



is given by E^ = 3/4E^ while Coulomb disorder scales 
as n^, where n is the number of electrons trapped in- 
side the QDs [2^. In order to fit the data we assume 
E^{B) — E^{0) and retained the temperature depen- 
dence (from Eq. [1]). Although an approximation, it is 
instructive to present this analysis as it represents an 
upper bound on the g-tensor, and most importantly, it 
can help in elucidating the symmetries for this particu- 
lar state. This reasonably simple model which takes into 
account the most important factors in determining the 
Zeeman splitting in quantum dots, including tempera- 
ture effects and wavefunction compression, describes well 
the whole data set. 

From the g factor obtained at each angle and 
for each shell, the g tensor modulus was de- 
termined for the polar and azimuthal scans by 
^^ cos(6')2 -H g2^^Qj sin(6l)2 and \gaz\ = 

sin(0)2. Figures 2 and 3 show the 
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polar and azimuthal scans. On the top panel the exper- 
imental set-up is represented showing the orientation of 
the magnetic field with respect to the QD crystalline axis, 
as well as the FD wavefunctions for each shell calculated 
for B = 15T. The data are shown on the bottom panel 
with the corresponding fits. 

From the polar scan one can immediately note that for 
the s shell the g-factor is quite anisotropic whereas for the 
p shell it is constant and always smaller. This is a quite 
surprising result at first, considering the the single par- 
ticle energy (htoo — 37.8meV) and confinement along z 
{Ez = — 280meV), and the fact that the z-component of 
the wavefunction for both states is basically the same. 
From the fit, we find Ifffgoiil = 1-51 ± 0.03 and |gM;^oll ~ 
0.57 ±0.05, and {g^^^^^l « {g^^^^^l = 0.47 ± 0.07. The 
larger values for g^ are consistent with a stronger local- 
ization inside the QD due to the deeper confinement. In 
order to draw a qualitative picture for the wavefunctions, 
we plot the solutions for the s and p wavefunctions for 
different magnetic field tilt angles (figures 2 and 3, top 
panel). Insofar as the behavior for Ig"] is concerned, one 
notes that the wavefunction is more localized into the 
QD and consequently more sensitive to the confinement 
potential details at the QD center. It is expected a high 
anisotropy given the pancake geometry of the QDs. The 
somewhat unexpected result comes about for the p shell. 
A highly isotropic behavior is found, and we associate 
this behavior to the symmetry of the p- wavefunction, 
as depicted in the upper panel of figure 2. First, as the 
electrons on the p shell are more delocalized, a leakage of 
the wavefunction along the growth direction takes place, 
which brings the g-factor modulus values closer to that 
of the matrix (IgoaAs — 0.44|). Second, the wavefunc- 
tion for the p shell has a node at the center of the QDs, 
which permits probing of the regions outside or at the 
interfaces of the QDs. Hence, by measuring the g-factor 
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FIG. 2: (top) Wavefunction for the s and p shell for three 
directions on applied magnetic field ^ = 0, ^ = 7r/3 and 6 = 
7r/2. x,y and z correspond to [110], [110] and [001] directions; 
a QD is schematically represented in the same plot, (bottom) 
g-tensor for the s and p shells for the polar scan. 



for the s and p shell, one can evaluate the details of the 
confining potential at selected spatial locations, in a sim- 
ilar fashion that was carried out for the determination of 
the chemical environment of deep levels [7| . 

Figure 3 shows the g-tensor for an in-plane field con- 
figuration (azimuthal scan). For both s and p shells the 
g-factor is independent of the field direction within the 
experimental uncertainties, consistent with a cylindrical 
symmetry for the QDs. 

Several competing effects have to be taken into account 
when interpreting the obtained results: i) local strain and 
local crystal fields which change the details of the confin- 
ing potential, ii) non-uniform composition, iii) quantum 
confinement, and iv) non-parabolicity, which takes place 
in InAs. All these parameters may influence the spin- 
orbit coupling, which is one of the important components 
in the g-factor determination. Modifications into g-factor 
has been demonstrated by tuning both strain 20|,|23| and 
composition [2j] in the QDs. If one compares these ex- 
perimental results with theory, a good agreement is found 
which corroborates the description by Pryor et al [21,] . In 
essence, we find that effects i-iii are basically the same, 
and they by and large determine the behavior for the g- 
tensor. More careful experiments on different samples are 
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FIG. 3: (top) Wavefunction for the s and p shell for three 
directions on applied magnetic field <^ = 0, = 7r/4 and (j) ~ 
■K 12. x,y and z correspond to [110], [110] and [001] directions; 
a QD is schematically represented in the same plot, (bottom) 
g-tensor for the s and p shells for the azimuthal scan. 



required to draw a more complete picture on the relative 
effect of each component. 
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FIG. 4: (left) calculated and measured addition energies for 
the polar scan; (right) same for the azimuthal scan. 

As a final verification of the models utilized in this 
work, we calculate the electron addition energy for n=2, 
i.e., s shell filhng, from the values obtained from the fits. 
The addition energy can be calculated according to equa- 
tion[TJ Figure [4] shows the calculated (left) and measured 



(right) results for the polar and azimuthal scans. The 
usefulness of this verification is primarily for evaluation of 
issues such as adequacy of the proposed model, signal to 
noise ratio, and an overall picture of what can be resolved 
for this two electron system. An important point men- 
tioned previously was that of ensemble measurements. In 
this experiment, one obtains a more reliable evaluation of 
the g-tensor as one is averaging over many different spin 
configurations which are temperature and magnetic-field 
dependent. This experiment thus provides a more repre- 
sentative description of g-tensor in nanostructures. 

In summary, we have inferred the g-tensor for the s 
and p shells of self-assembled QDs. We found that for 
the s shell the g-tensor is highly anisotropic, reflecting 
the confinement potential details. For the in-plane com- 
ponent and for the p shell, the modulus of the inferred 
g factors were close to the bulk GaAs value. Finally, 
we found that for the p shell the g tensor was isotropic 
within our experimental resolution, which is consistent 
with wavefunction having a node at the QD center and 
being more delocalized along the growth direction. We 
acknowledge the financial support by CNPq, FAPESP 
and HP Brazil, and the usage of the high magnetic field 
facihty at GPO/IFGW - UNICAMP. 
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